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Abstract 

Oscillating behaviour of the susceptibility x an d heat capacity C is considered 
for normal and superconducting mesoscopic systems (nanoclusters and quantum 
dots). It is proved that at low temperature an increasing magnetic field applied to 
a mesoscopic system generates local extrema of x an d C. A maximum for x an d a 
minimum for C simultaneously arise in those points of the field where crossings of 
quantum levels of the normal and superconducting mesoscopic systems take place. 
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1 Introduction 



Low temperature oscillations of diverse electromagnetic and thermal charac- 
teristics (resistivity, magnetization, susceptibility, thermal conduction, heat 
capacity etc.) of macroscopic and mesoscopic bodies in an increasing uniform 
magnetic field are apparently a general property of systems formed by charged 
particles. These oscillations have been last years studied experimentally and 
theoretically in mesoscopic systems of different geometry and microscopical 
structure: in mesoscopic rings [1], in a thin spherical layer [2], in superconduct- 
ing discs [3,4,5,6,7], in micron-size heterostrucures [8]. Observing such variety 
of properties and objects one is forced to arrive at a conclusion that there must 
be a common reason for arising this phenomenon. In this Letter as a source 
of such oscillations we regard crossings of the quantum levels of a mesoscopic 
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system abstracting from the physical nature of the levels and system itself. 
We concentrate here on static properties (the magnetic susceptibility x an d 
heat capacity C) in the regime of the thermodynamical equilibrium of the 
canonical ensemble at low temperatures. 

We assume that a system exposed to a uniform magnetic field (B) possesses 
a number of discrete quantum states altering their mutual energy disposition 
depending on B. Thereby for each value of B there exist a lowest state, the 
ground state for a given B, and the nearest to it first excited state with 
energies E (B) and Ei(B) respectively. Each of these energies is a smooth 
function of B, however E Q (B) and Ei(B) can locally have opposite slopes 
with B. Therefore the difference E W (B) 

E W (B) = E 1 (B) - E (B); E w (B o ) = 0. (1) 
can become equal to zero for a value of B = B . 

In Sec. 2 we show that this point of the level crossing is a peculiar point where 
at very small temperatures both x and C reach local extrema: a maximum for 
X and a minimum for C. 



2 Extrema in x an d C 

In the low temperature limit only two first terms practically exhaust the canon- 
ical partition function 

Z = exp[-pE (B)} + exp[-pE l (B)} + Z v , (3 = 1/T. (2) 

The rest, Z T , gives a weak background for x an d C near the point B in the 
vicinity of which the partition function can be expressed via a small parameter 

Z = ±{exp[-PE w (B)]-iy, = 0, (3) 

\nZ = -(3E (B)+\n(l + Z) + a + \n2, 
a = ]-J2 ex P{-P(E l -E )}<^l. 

1 i>2 

The magnetic susceptibility and heat capacity are naturally divided near B 
into smoothly varying parts Xs, C s and the rapidly varying functions Xe, C e 
with extrema in Bq\ 
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C(B )= C s + C e , (6) 
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V is the volume of the system. The derivative 8/ dBo implies hereafter the 
derivative with respect to B in the point of B — B . 

At low temperatures (T < 5, 5 is the mean level spacing at B = 0) we take 
into account only those terms in Eqs.(5) and (7) which are proportional to 
the maximum power of (3 at B = B , <9™£/&B™ ~ 0.5(-f3) n (dE w /dB ) n , that 
gives 
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Xc is paramagnetic (% s determines the diamagnetic component of x) an d pro- 
portional to the squared difference of the magnetic moments of the crossing 
states. Thus, Eqs.(8) -(11) confirm the existence of a local maximum in x 
and a minimum in C at the level crossing point B . Increasing the field can 
generate appearance of new level crossing points that will give x an d C an os- 
cillating functions of B. Intervals between peaks and their amplitudes depends 
on the physical nature of the system and its shape. Here we limit ourselves 
by considering weak fields (the cyclotron radius is much larger than the size 
of the system). The alterations of quantum levels v.s. B are caused by remov- 
ing the spin and rotational degeneracy, the latter occurs when a system has 
rotational symmetry. To investigate oscillating behavior of x an d C in Sec. 3 
we apply the independent electron model to normal mesoscopic systems and 
in Sec. 4 superconducting clusters of the perfect spherical shape are described 
by using the exact solution of the conventional superconducting hamiltonian. 
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3 Magnetic oscillations in normal clusters 



The character of oscillations in electron \ and C of a mesoscopic system is 
dependent on the its shape. In spherical clusters at B — each many-fold 
degenerate level (the spin-orbital interaction is omitted here) displays a sheaf- 
wise splitting. Hence magnetic sublevels originated from one spherical /-shell (/ 
is an orbital momentum) cannot cross each other with the growth of B. Thus 
the minimum strength of a magnetic field, B min , (here we consider the case 
of the closed Fermi-shell at B = 0) required to arouse the first crossing and 
the first extremum is determined by 5 sn , the distance between the Fermi-shell, 
Zi^-shell, and the unoccupied Z^+i-shell at B = 0: fi* B B m i n (lF + If+i + 2) = <5 sn , 
fi* B = ji B m* /m, m* is the effective electron mass. We assume that orbital and 
spin projections are aligned at the first crossing. Nevertheless the absolute 
value of -Bmin is practically independent of /, since the level splitting near 
F at B = is proportional to I: 5 sn ~ (21 p + 1)5, where 5 = Asf/^N is 
the mean level spacing of the Fermi gas near the Fermi energy (s.f). Thus, 
-Bmin ~ Sm*/fiB m and for nanometer normal metallic clusters (N ~ 10 5 , 
Bp ~ 10 eV and m — m*) B min is of the order of IT. The evolution of some 
single-electron levels in a spherical cavity and corresponding extrema in x and 
C are shown in Fig.l. 

Properties of quantum dots are frequently interpreted in the two-dimensional 
cylindrically symmetric oscillator model [9] . This symmetry also gives a rather 
high degeneracy of electron levels near the Fermi level, which is removed in 
the same manner as in the previous spherical case. B m \ n of the first level 
crossing is assessed analogously: B m \ n ~ e^m* /2iV /jgm, N is the free electron 
number in the quantum dot. As Ep and m*/m for such systems can be small 
(e.g. for GaAs m*/m ~ 0.07, Bp ~ 10 meV [9]) oscillations of x could be 
observable in relatively weak fields for a wide range of N. x and C v.s. the 
cyclotron frequency are shown in Fig. 2, 3 for N ~ 10 3 ( we consider the case 
of the closed cylindrical Fermi-shell at T = B = ). The feature of x and C 
oscillations in systems with cylindrical or spherical symmetry is the presence 
of several groups of peaks corresponding firstly crossings of levels from two 
nearest cylindrical or spherical shells separated at B = by one shell gap, 
then cross levels separated at B = by two shell gaps and so on. Increasing 
temperature (Fig. 3) leads to confluence of separate peaks so that curves for x 
and C become similar and their gross structure is caused by the shell structure 
of the single-particle energy spectrum. 

Much simpler behavior of x and C as function of B is observed in clusters 
having no rotation symmetry axis. In this case level crossings are caused by the 
Zeeman splitting and perturbation of orbital motion by weak fields proves to 
be insignificant, because orbital energy shifts for such systems are quadratic in 
B. Then signs of energy denominators appearing in perturbative calculations 
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are opposite for levels near the Fermi-level that reduces second order terms. 
The example of such "spin" oscillations of x is concerned with electron system 
in the oscillator well with oscillator frequencies Q x = 0.9Q y = 0.5Q Z that 
corresponds to an ellipsoid with semiaxes a x ~ l.la y = 2a z (we adopt the 
relation between semiaxes and oscillator frequencies from Ref. [10]: fl x a x = 
Qydy = VL z a z ). Such system has no rotation axis and a magnetic field directed 
alone; generates "spin" oscillations (Fig. 4). These oscillations of x i n 

the anysotropic oscillator are much like oscillations in the equal level spacing 
model [11]. 



4 Magnetic oscillations in the gapless superconducting region 

Another example of magnetic oscillations of x an d C is given by the single- 
shell model (SSM) of superconductivity. Such model can be a reasonable ap- 
proach to a description of spherical superconducting clusters with numbers 
of conduction electrons 10 3 < N < 10 5 [12]. There is a wide range of orbital 
momentum lp for states successively becoming the Fermi level with growth 
of N. Averaging through N gives l sh ~ A 1 / 3 [12]. At B = the level spacing 
between degenerate shells 5 s h ~ 21 F 5 ~ 8e>/3 A 2 / 3 is much more than the 
BCS superconducting gap at B — T — 0, A B cs(0), e.g. 5 sh > 10 meV and 
Abcs(0) < 1 meV for AZ-clusters with iV under consideration. Therefore the 
superconducting pair correlations can arise between electrons (their number 
is iV s h) occupying only one shell with a large /. 

The Hamiltonian of SSM consists of the pairing interaction with the strength 
G and the interaction with a field B. The exact spectrum of such Hamiltonian 
is well known [10] and at B = is determined by the seniority v (unpaired 
electron number) 

E v = -j{N sll - u)(Al + 4 - A sh - u). (12) 

For each v there is a multiplet of states (degenerate at B — 0) with different 
orbital and spin projections which are splitted by a field. The energy of the 
lowest state for each v is 

E„(u>) = E V - + 2 - u) - \[l - R1 + 4(1 - 5 U . )}, (13) 

Thus with growth of B (in this Section uo = hujfi* B B) the role of the ground 
state is successively played by states with v = 0,2,4,... for even N sh or 
v = 1,3,5. . . for odd A sh , i.e. increasing B leads to a sequence of crossings 
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of the i/- level with v + 2-level, resulting in a step-wise decrease of the pairing 
gap (the region of the gapless superconductivity). 

To take into account the temperature dependence of x an d C quantities we 
have analytically summed up those components of the canonical partition 
function Z which involve states inside the considered /-shell 



Z = J2exp(-f3E u )[<S> u - $„_ 2 (1 - 8„. )], (14) 



u/2 2 



rT„:.. 1 + <V.O 



Vmin = |[1 - (-)"], 



£ sinh(/^§) ' ^ k B r 



Fig. 5 shows x an d C as functions of the magnetic field and temperature. The 
latter is measured in T cr i t , the critical temperature of the conventional BCS 
theory for half-filled shell. Oscillations of x an d C created by level crossings in 
SSM spread on a range of rather small fields u> ~ Abcs(0)//f [12] (B < 0.5 T 
for AZ-clusters with iV > 10 3 and m — m*). In nonspherical superconducting 
clusters where a weak magnetic field does not influence orbital motion, the 
level crossings are determined only by alignment of spins. Therefore greater 
values of B min are expected to induce the first extrema in x an d C. 



5 Conclusion 



In conclusion we have shown that an increasing uniform magnetic field applied 
to a mesoscopic system transforms the structure of the ground state in a 
such a way that in level crossing points synchronously appear local extrema 
in magnetic susceptibility (maxima) and heat capacity (minima). Thereby 
these thermodynamic quantities oscillate with growth field. Periodicity and 
amplitudes of the low temperature oscillations carry information concerning 
properties of the quantum states of investigated normal and superconducting 
systems. 
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hoo /8 

c 

Fig. 1. Electron levels vs the cyclotron frequency (fiuj c = 2/j,* B B, fi* B = fiB m * / m ) 
and magnetic oscillations of the susceptibility \ (dotted line) and heat capacity C 
(solid line) for a sphere with TV" ~ 10 5 . Points of level crossings where x anci C reach 
local extrema are marked by solid circles. 5 = 4ef/3N. The bold line corresponds 
to the Fermi-level. \ is given in units 10 • | xl I, Xl = —^%N/2Vsf is the Landau 
diamagnetic susceptibility of the degenerate free electron gas, m = m*. 
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Fig. 2. Low temperature magnetic oscillations \ (dotted line) and C (bold line) in 
a two-dimensional cylindrical oscillator with N ~ 10 3 . 8 = Aef/3N. \ ls given in 

XL- 
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T/5 = 0.2 




I i I i I i I i I i I i L 

1 2 3 4 5 6 

h co/5 



Fig. 3. x vs the cyclotron frequency to c at different temperatures for a 
two-dimensional cylindrical oscillator with N ~ 10 3 (closed Fermi-shell) 
5 = Ae F /3N. 
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Fig. 4. x vs w c at different temperatures for an anysotropic oscillator 
(Q x = 0.9Q y = 0.5O 2 ) with N ~ 10 5 , magnetic field is directed along x -axis. 
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Fig. 5. Magnetic oscillations in a superconducting sphere. B cr [ t is the critical mag- 
netic field of the BCS theory at T = 0, A BCS (T = 0, B crit ) = 0. (a) x vs B (xd is 
the diamagnetic susceptibility for an ^/-cluster with N = 10 3 ). (b) C vs B for the 
same cluster. 
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